Large-scale magnetic topologies of M dwarfs 
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Abstract. We present here the first results of a spectropolarimetric analysis of a small sample (~ 20) 
of active stars ranging from spectral type MO to M8, which are either fully-convective or possess 
a very small radiative core. This study aims at providing new constraints on dynamo processes in 
fully-convective stars. 

Results for stars with spectral types M0-M4 - i.e. with masses above or just below the full 
convection threshold (^0.35 Mq) - are presented. Tomographic imaging techniques allow us to 
reconstruct the surface magnetic topologies from the rotationally modulated time- series of circularly 
polarised profiles. 

We find strong differences between partly and fully convective stars concerning magnetic field 
topology and characteristic scales, and differential rotation. Our results suggest that magnetic field 
generation in fully convective stars relies on different dynamo processes than those acting in the 
Sun and other partly convective stars, in agreement with theoretical expectations. 
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CONTEXT 

In partly convective stars such as the Sun, magnetic fields - the energy source of most 
activity phenomena - are induced by plasma motions: the combined action of differential 
rotation (Q. effect) and cyclonic convection (a effect) manages to generate a self- 
sustained magnetic field. The so-called aQ. dynamo [ ] is believed to operate mostly 
through the tachocline, a thin zone of strong shear located at the interface between the 
radiative inner zone and the convective envelope [ ]. 

Stars with masses lower than about 0.35 Mq are fully-convective [ ] and thus do 
not possess a tachocline. However, they manage to trigger magnetic fields [4, 5, 6] 
and are very active [7, 8, 9]. Though significant progress was made since first non- 
solar dynamo mechanisms were proposed [ ], theoretical and numerical modelling 
require observational constraints. It is now acknowledged, from observational [11, 12] 
and theoretical points of view [13, 14, 15], that fully convective stars (FCS) manage to 
yield large scale magnetic fields. But the properties of such magnetic fields, and their 
dependency on stellar parameters (in particular mass and rotation rate) are not yet clear. 

We present here the first results of a spectropolarimetric analysis of a small sample 
of active M dwarfs with spectral types ranging from MO to M8, which are either fully 
convective or possess a very small radiative core. We aim at exploring the properties 
of the large-scale magnetic topologies of FCS, and their evolution with main stellar 
parameters. The stars were selected from the rotation-activity study [ /]. We chose only 



TABLE 1. Fundamental parameters of the stellar sample. Columns 1-4 respectively list 
the name, the spectral type, the stellar mass, and the logarithmic relative X-ray luminos- 
ity logRx =log(Lx/Liyoi). The projected rotation velocity and rotation period inferred from 
Zeeman Doppler Imaging (ZDI) are mentioned in columns 5 and 6. Columns 7-10 respec- 
tively list the empirical convective turnover time from [ ], the effective Rossby number 
p 

Ro = the theoretical radius suited to the stellar mass from [ ], and the inclination 
angle used for ZDI. See [18, 19] for more details. 
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active stars so that the magnetic field is strong enough to produce detectable circularly 
polarised signatures, allowing us to apply tomographic imaging techniques. More details 
about the sample are available in Table 1 . 

For this study, we used the twin instruments ESPaDOnS on the 3.6-m Canada- 
France-Hawaii Telescope (CFHT) located in Hawaii and NARVAL on the 2-m Telescope 
Bernard Lyot (TBL) in southern France. These spectropolarimeters, built on the same 
design, can produce Stokes I,Q,U and V spectra spanning the entire optical domain 
(from 370 to 1000 nm) at a resolving power of ^ 65000 [ ]. 

We performed monitoring observations of the sample in circularly polarised light 
(Stokes V). Least-squares deconvolution (LSD) [ ] was then applied, resulting, for 
each spectra, in a synthetic line profile gathering polarimetric information from most 
photospheric atomic lines. 

IMAGING PROCEDURE AND MODEL DESCRIPTION 

For each star of the sample, our aim is to infer the topology of the surface magnetic 
field from the circularly polarised (Stokes V) LSD profiles we obtained. This can be 
achieved using a Zeeman-Doppler Imaging (ZDI) code [22]. The imaging process is 
based on the principles of maximum entropy image reconstruction. The magnetic field 
is decomposed into its poloidal and toroidal components, both expressed as spherical 
harmonics expansions. Starting from a null magnetic field, we iteratively improve our 
magnetic model by comparing the synthetic Stokes V profiles with the observed LSD 
profiles, until we reach an optimal field topology that reproduces the data at a given 
level. The inversion problem being partly ill-posed, we use the entropy function to select 



the magnetic field with lowest information content among all those reproducing the data 
equally well [ ]. 

To compute the synthetic Stokes V profiles, the star is divided into a grid of ^ 
1000 cells on which the magnetic field components are computed directly from the 
coefficients of the spherical harmonics expansion. The contribution of each individual 
pixel is computed from a model based on Unno-Rachkovsky's equations [ ]. We then 
integrate all contributions from the visible hemisphere at each observed rotation phase. 

While computing the Stokes V profiles it is possible to account for differential rota- 
tion. For a given differential rotation law, each local line profile is Doppler- shifted as a 
function of the observation phase. It is then possible to investigate how the fit quality 
varies in a reasonable range of Prot and dQ. values. We can thus derive the optimal Prot. 
dQ. and corresponding error bars [24, 25, 12]. 

RECONSTRUCTED MAGNETIC FIELDS 

For each star of our sample, we reconstruct with ZDI the large-scale surface magnetic 
field at least up to degree £ = 6. We also measure differential rotation and assess time- 
variability of the magnetic topology whenever possible. 

The field is characterised by three quantities: (a) the overall magnetic energy, (b) 
the ratio of magnetic energy reconstructed in the poloidal modes and (c) the ratio of 
poloidal magnetic energy reconstructed in the axisymmetric modes (defined by m < |). 
The results are described below and presented in a more visual way in Figure 1 as a 
function of and Prot- See [Donati, these proceedings] for a version of this diagram 
including G and K dwarfs. To compare magnetic field generation in stars of different 

masses, it is convenient to introduce the effective Rossby number which rescales Prot by 

p 

a mass dependant coefficient. Ro = where is an empirical convective turnover 
time inferred from X-ray luminosities [16]. On Fig. 1 we also plot contours of constant 
Rossby number Ro = 0.1 and 0.01 respectively corresponding approximately to the 
saturation and super-saturation thresholds [e.g., 26]. 

The magnetic field topologies we reconstruct with ZDI and the differential rotation 
amplitudes are very different on each part of the O.SMq boundary. This threshold 
is very sharp and well defined, with little apparent dependence with the rotation period. 
Typical examples of magnetic maps on each part of this limit are shown in Fig. 2. 

For stars more massive than 0.5 M©, we recover magnetic topologies including 
(i) a strong toroidal component and (ii) a high non-axisymmetric degree of the 
poloidal component. For 4 stars of this subsample we can derive differential rotation. 
We find dQ. > dQ.Q, in agreement with the dispersion observed in previous photometric 
measurements of Prot- As a consequence, surface magnetic features are short-lived, 
the topology completely changes on a timescale of a few months. 

In the low mass subsample (M^<0.5 M©), we reconstruct much stronger magnetic 
fluxes and very different large-scale magnetic topologies : (i) mostly poloidal (~ 90% 
of the reconstructed energy), (ii) strongly axisymmetric (except for EV Lac, more than 
half of the magnetic energy is reconstructed in m = modes) and (iii) close to a dipole 
with more 50% of the reconstructed magnetic energy lying in poloidal modes of degree 




Prot (d) 

FIGURE 1. Properties of the magnetic topologies of M dwarfs as a function of rotation period and 
stellar mass. Larger symbols indicate larger magnetic fields while symbol shapes depict the different 
degrees of axisymmetry of the reconstructed magnetic field (from decagons for purely axisymmetric fields 
to sharp stars for purely non axisymmetric fields). Colours illustrate the field configuration (dark blue 
for purely toroidal fields, dark red for purely poloidal fields and intermediate colours for intermediate 
configurations). Solid lines represent contours of constant Rossby number Ro = OA and 0.01. The 
theoretical full-convection limit (M^ O.SSM© [ ]) is plotted as a horizontal dashed line. 



£ = 1 . Very weak differential rotation is inferred from our data with 3 stars having dQ. 
of the order of a few mrad d~^ This very weak differential rotation is in agreement with 
the most recent numerical simulations These stars were observed at two different 
epochs separated by ~ 1 yr. Evolution of the magnetic topologies is small, in some 
cases, it is possible to fit observations separated by ~ 1 yr with a unique magnetic 
topology. 

CONCLUSIONS AND PERSPECTIVES 

Very different large-scale magnetic topologies are observed on each part of the 
0.5 M0 limit. We also note that dynamo processes become suddenly much more efficient 
at triggering large-scale magnetic fields (see Fig. 1 and 3) at approximately the same 
mass (c^ 0.4 Mq). This strong step is not visible in the log7?x vs Ro plot. This result 
suggests that (i) the X-ray emission is sensitive to overall magnetic energy whereas we 
are only sensitive to the largest scales, (ii) At a given Ro stars with mass above or below 
^0.5 M0 generate comparable magnetic energy, but with very different spatial scales 
repartition (the less massive stars triggering more magnetic energy in the largest scales). 
These strong changes in magnetic field generation occur at masses slightly larger 
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FIGURE 2. Upper row: surface magnetic flux of DT Vir (0.59 M©) as derived from our 2008 data set. 
The three components of the field in spherical coordinates are displayed from left to right (flux values 
labelled in G). The star is shown in flattened polar projection down to latitudes of —30°, with the equator 
depicted as a bold circle and parallels as dashed circles. Radial ticks around each plot indicate phases of 
observations. Bottom row: same for YZ CMi (0.31 Mq). 
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FIGURE 3. Left panel: Reconstructed magnetic energy as a function of Ro. Right panel: logarithmic 
relative X-ray luminosity as a function of Ro (see Tab. 1). 



than the theoretical hmit to full convection {c^ 0.35 M©). This may be due to strong 
shrinking of the radiative inner zone predicted by theoretical models, from ^ 0.5 at 
= 0.5 M0 to nothing at 0.35 Mq [17, 27]. 
We are currently completing this survey to extend it to fast rotators with > 0.5 M©, 
slow rotators with 0.2 <M^< 0.5 M©, two regimes not explored in the present sample. 
Present efforts are also directed to ultracool dwarfs (M^ < 0.2 M©), to investigate how 
dynamo processes operate down to the brown dwarf limit, i.e. when stellar atmospheres 
start to become neutral. 
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